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Abstract 

In this paper, we investigate optimal control problems for Allen-Cahn equations with 
singular nonlinearities and a dynamic boundary condition involving singular non- 
linearities and the Laplace-Beltrami operator. The approach covers both the cases 
of distributed controls and of boundary controls. The cost functional is of standard 
tracking type, and box constraints for the controls are prescribed. Parabolic prob- 
lems with nonlinear dynamic boundary conditions involving the Laplace-Beltrami 
operation have recently drawn increasing attention due to their importance in ap- 
plications, while their optimal control was apparently never studied before. In this 
paper, we first extend known well-posedness and regularity results for the state equa- 
tion and then show the existence of optimal controls and that the control-to-state 
mapping is twice continuously Frechet differentiable between appropriate function 
spaces. Based on these results, we establish the first-order necessary optimality 
conditions in terms of a variational inequality and the adjoint state equation, and 
we prove second-order sufficient optimality conditions. 



1 Introduction 

Let Q C JR N , 2 < iV < 3 , denote some open and bounded domain with smooth boundary 
T and outward unit normal n , and let T > be a given final time. We put Q := Q x (0, T) 
and S := T x (0,T) , and we assume that Pi > 0, 1 < % < 6, are given constants which 
do not all vanish. Moreover, we assume: 
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(Al) There are given functions 

z Q eL 2 (Q), z E GL 2 (E), z T eH\Q), zrpeH^r), 
Ux,U2 G L°°(Q) with u\ < U2 a. e. in Q, 
Ui r ,u 2r G L°°(£) with U\ T < u 2r a.e. m ^ • 

We then consider the following (tracking type) optimal control problem: 
(CP) Minimize 

J((v,Vr),(u,u T )) := y / f \y - z Q \ 2 dx dt + ^ [ f \Vr ~ z^\ 2 dV dt 

+ ^Jjy(.,T)-Z T \ 2 dx + P±j\y F (.,T)-Z r>T \ 2 dT 

+ T /7 l M l 2dxdt + V l^l 2 drdt (1.1) 
^ Jo Jn ^ y y r 

subject to the parabolic initial-boundary value problem with nonlinear dynamic boundary 
condition 

y t -Ay + f'(y)=u a.e. in Q, (1.2) 
<9i?/r - A r yr + <9 n y + ^'(yr) = «r, y\r = Vv, a. e. in S, (1.3) 
2/('i0) = y o a.e. in Q, y r (-,0) = ?/or a.e. on T, (1.4) 
and to the control constraints 

(u, ur) G W a d := G L 2 (Q) x L 2 (S) : 

ui < w <U2 a.e. in Q, «i r < wr < «2 r a.e. in E } . (1.5) 

Here, yo and yo r are given initial data, Ap is the Laplace-Beltrami operator on T , and 
the functions / , g are given nonlinearities, while u , up play the roles of distributed or 
boundary controls, respectively. Note that we do not require Mr to be somehow the 
restriction of a on T; such a requirement would be much too restrictive for a control to 
satisfy. 

We remark at this place that for the cost functional to be meaningful it would suffice to 
only assume that zt G L 2 (Q) and z^.t G L 2 (T) . However, the higher regularity of zt 
and z Tt requested in (Al) will later be essential to be able to treat the adjoint state 
problem. 

The system fll.2p - fll.4p is an initial-boundary value problem with nonlinear dynamic 
boundary condition for an Allen-Cahn equation. In this connection, the unknown y 
usually stands for the order parameter of an isothermal phase transition, typically the 
fraction of one of the involved phases. In such a situation it is physically meaningful to 
require y to attain values in the interval [0, 1] on both Q and T . A standard technique to 
meet this requirement is to postulate that the first derivatives of the bulk potential / and 
of the surface potential g become singular at and at 1. A typical form for such a poten- 
tial is / = fi + /2 , where f 2 is smooth on [0, 1] and fi(y) = a [yln(y) + (l — y) ln(l — y)] 
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with some a > . Another possibility is to choose f± as the indicator function i^i] of 
the interval [0, 1] ; in this case f[ has to be replaced by the subdifferential <9/[ 0i i] , and 
/ becomes a double obstacle potential. In this case, (11. 2p has to be understood as a dif- 
ferential inclusion or variational inequality. Similar choices can be made for the surface 
potential g. 

There exists a vast literature on the well-posedness and asymptotic behaviour of the Allen- 
Cahn equation with the no- flux boundary condition d n y = in place of (11.31) . Also, the 
well-posedness and asymptotic behavior of the system (ll.2l) - (ll.4l) has been the subject of 
numerous papers (see [3] and the many references given there). 

Moreover, distributed and boundary control problems for the Allen-Cahn equation with 
no-flux boundary conditions or boundary conditions of the third kind have been studied 
in a number of recent papers, in particular, for the case of the double obstacle potential. 
In this connection, we refer to j5] and [6]. Associated stationary, that is, elliptic MPEC 
problems have been studied in [5] (see also the monograph [2]), and the related Cahn- 
Hilliard case was recently analyzed in (TUJ. We also like mention the works [7], [8], [12] 
that treated optimal control problems for the Caginalp-type temperature-dependent gen- 
eralization of the Allen-Cahn equation in the case of nonsingular potentials and standard 
boundary conditions; a thermodynamically consistent temperature-dependent model with 
singular potential of the above logarithmic type was the subject of [IB] . 
The main novelty of the present paper is to study optimal control problems with singular 
potentials of the logarithmic type and dynamic boundary conditions of the form (11.31) . 
In fact, while various types of dynamic boundary conditions have already been studied 
in connection with optimal control theory (see [IT] , for a recent example), it seems that 
dynamic boundary conditions involving the Laplace-Beltrami operator have not been 
considered before. One of the difficulties is that from the viewpoint of optimal control it 
does not make sense to postulate that the controls u and u-p satisfy U|r = «r • 
The paper is organized as follows: in Section 2, we give a precise statement of the problem 
under investigation, and we derive some results concerning the state system (ll.2l) - (ll.4p 
and a certain linear counterpart, which will be employed repeatedly in the later analysis. 
In Section 3, we then treat the optimal control problem, proving the existence of optimal 
controls and deriving the first-order necessary and the second-order sufficient optimality 
conditions. During the course of this analysis, we will make repeated use of the elementary 
Young's inequality 

ab < 7|a| 2 + — \b\ 2 Va,6e H V7 > 0, 
47 

of Holder's inequality, and of the fact that we have the continuous embeddings if x (f2) C 
L P (Q) , for 1 < p < 6, and H 2 (Q) C L°°(Q) in three dimensions of space. In particular, 
we have 

\\v\\lp(q) < C p \\v\\ H i {n) Vt; e H\n), (1.6) 
IMU<»(n) < C*oo |M|fl3(n) VveH 2 (tt), (1.7) 

with positive constants C p , 1 < p < 00 , that only depend on Q . 
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2 General assumptions and the state equation 

In this section, we formulate the general assumptions of the paper, and we state some 
results for the state system fll.2p - fll.4l) . To this end, we introduce the function spaces 

H:=L 2 (n), V-.^H 1 ^), H r :=L 2 (T), V r := H^T), 

H := L 2 {Q) x L 2 (S), X := L°°(Q) x L°°(£), 

y-.= {(y,y r ) : y G i/ 1 (0, T; H) fl C°([0, T];V)C\ L 2 (0, T; H 2 (Q)), 

y r G H\0, T; H r ) fl C°([0, T]; Vb) n L 2 (0, T; 2 (r)), y r = y\r} , (2.1) 

which are Banach spaces when endowed with their natural norms. In the following, we 
denote the norm in a Banach space E by || ■ \\e] for convenience, the norm of the space 
H x H x H will also be denoted by || • \\h ■ Identifying H with its dual space H* , we 
have the Hilbert triplet V C H C V* , with dense and compact embeddings. Analogously, 
we obtain the triplet Vp C H r C Vp* , with dense and compact embeddings. We make the 
following general assumptions: 

(A2) / = /i + / 2 and g = g 1 + g 2 , where f 2 , g 2 G C 3 [0, 1] , and where f 1 ,g 1 G C 3 (0, 1) 
are convex and satisfy the following conditions: 

lim f[(r) = lim g[(r) = —00 , lim f[(r) = lim g[ (r) = +00 . (2-2) 

r\0 r\0 r/*l r /*1 

3 Mi > 0, M 2 > such that |/J(r)| < Mi + M 2 |^(r)| Vr G (0,1). (2.3) 

(A3) y eV , y 0r eV r , and we have /i(y ) G L x (fi) , g-i(y 0r ) ^ ^ : ( r ) , and 

< y < 1 a. e. in fi, < y 0r < 1 a - e. on T . (2.4) 

Remark 1: The condition (12.21) is obviously satisfied if f\ and g\ are potentials of 
logarithmic type as those mentioned in Section 1, while (12.31) is needed for the existence 
result from [3] that will be used below. 

To simplify notation, in the following we will denote the trace y\? (if it exists) of a function 
y on T by yr without further comment. Now observe that set U a d is a bounded subset 
of X . Hence, there exists a bounded open ball in X that contains W a( j • For later use it 
is convenient to fix such a ball once and for all, noting that any other such ball could be 
used instead. In this sense, the following assumption is rather a denotation: 

( A4) U is a nonempty open and bounded subset of X containing W at } , and the constant 
R > satisfies 

IMU°°(Q) + ll n r||L°°(s) < R V (u,u T ) eU. (2.5) 

The following result follows as a special case from [31 Theorems 2.3-2.5 and Remark 4.5] 
if one puts (in the notation of [3]) (3 = f{ , (3? — 9i > 71 — f% , n r = 92 there. 
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Theorem 2.1 Suppose that the general assumptions (A2), (A3) are satisfied. Then we 
have: 

(i) The state system (ll.2p — (jl.4p has for any pair (u, ur) El-La unique solution (y, y r ) E y 
such that 

< y < 1 a. e. in Q, < yo r < 1 a. e. on E . 

(ii) Suppose that also ( A4) is fulfilled. Then there is a positive constant K{ , which only 
depends on Q , T , y , y 0r , f , g, and R, such that for every (u,ur) the associated 
solution (y, yr) E y satisfies 

\\(y,yr)\\ 

+ \\yr\\H 1 (0,T;H r )nC ([0,T];V r )nL ,2 (0,T;H 2 (r)) < K{ , (2.6) 
\\f\y)\\L^T;H) + \W(yr)\\LH0,T;H r ) < K\. (2.7) 

Moreover, there is a positive constant K 2 , which only depends on Q , T , y , y 0r , f , 
g , and R, such that the following holds: whenever (u\,ui r ), (u2,^2 r ) E U are given and 
{yiiVi r )i (2/2, 2/2 r ) E y denote the associated solutions of the state system, then we have 

hi - y2\\c°([0,T];H) + H V (?/1 -y^Wl^Q) + llyir -2/2 r ||cO([0,T];if r ) 
+ ||V r (3/i r - 2/2 r )|||2 (S) 

< K* 2 hm - u 2 f L 2 {0>T . H) + \\u lr - u 2r \\ 2 L 2 (QiT . Hr) \ ■ (2.8) 



Remark 2: (i) It follows from Theorem 2.1, in particular, that the control-to-state 
mapping S , (u, tip) | — > S(u,ur) '■= (y,yr) is we ll defined as a mapping from X into 3^; 
moreover, S is Lipschitz continuous when viewed as a mapping from the subset U of % 
into the space (C°([0, T}; H) f) L 2 (0, T; V)) x (C°([0,T}; H r ) f] L 2 (0,T;V r )). 

(ii) Observe that we cannot expect y to be continuous in Q , since both d t yr and 
A r y r only belong to L 2 (S), so that also only d n yr E L 2 (S). However, we have y E 
L 2 (0,T;C°(Tl)) and y T E L 2 (0,T;C°(T)) . 

The next result is concerned with a linear problem with dynamic boundary condition. It 
will later be needed to ensure the solvability of a number of linearized systems. 

Theorem 2.2 Suppose that functions (u,ur) E H , C\ E L°°(Q) , c% E L°°(S) ; wq E 
if 1 (f2), and Wo r E H 1 ^) are given. Then we have: 
(i) The linear initial-boundary value problem 

w t — Aw + ci (x, t) w = u a. e. in Q, (2.9) 

d n w + d t wr — A r wr + c 2 (x,t) wr = u r a. e. on S, (2-10) 
w( ■ , 0) = Wo a. e. in Q, wr( ■ , 0) = wo r a. e. on T, (2-H) 
has a unique solution (w,Wr) E y . 
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(ii) There exists a constant C > 0, which only depends on Q, T, \\C\\\l°°(q) , and 
||c2||l°°(s) ; such the the following holds: whenever wq = and wo r = then 

\\(w,w r )\\ y < d\\(u,u r )\\ n . (2.12) 

Proof: We put j3{w) := (3r(w) := w and define the operators 

U(w)(x, t) := Ci(x, t)w(x, t) — w(x, t), H r (wr)(x, t) := c 2 (x, t)wr(x, t) — wr(x, t). 

With these definitions, we may rewrite the equations (12. 9p and (I2.10p in the form 

w t - Aw + P(w) + n(tw) = u, (2.13) 

d n w + d t w r - A r u>r + Pr{wr) + n r (u>r) = «r, ( 2 -14) 

respectively. Since the functions /3 and /?r are strictly monotone increasing in H, the 
system (12.111) . (I2.13p . (I2.14p has almost the same form as the system considered in Theo- 
rem 2.5 in [3], the only difference being that here II, Hr are linear and continuous operators 
while 7r, 7Tr in [3] were Lipschitz continuous functions. However, a closer inspection of 
the proof of Theorem 2.5 in [3] reveals that the argumentation used there carries over 
to the present situation with only minor and obvious modifications. Hence, the asserted 
existence result is valid. 

Now let Wq = and Wq f = . In the following, we denote by C{ , i G IN , positive 
constants that only depend on the quantities mentioned in the assertion of (ii). Testing 
(I2.13P by w t yields for every t G (0, T] the inequality 



f f iv 2 t dxdt + l - \\w(t)\\ 2 v + f f |^ r | 2 drdt + l - \\w T (t)\\ 



2 

Vr 



< I /((|ci| + l)M + |u|)K|drdt + / /((|c 2 | + l)|w r | + \u v \)\d t w v \&T &t , 
Jo Jn Jo Jr 

whence, using Young's inequality and the fact that c\ G L°°{Q) and c 2 G L°°(S), we 
obtain 

t r ft 

2a„A+ , IL..AAII2 , / / |Q„„ \2 -s-r A . , II.. /.mi 2 



wfdxdt + ||w(t)||v + / / \d t w r \ dTdt + \\w r 



t r rt 



<Ci(y j (\w\ 2 + \u\ 2 )dxdt + J J (\w v \ 2 + M 2 )drdt) . 

Gronwall's lemma then yields that 

\\w\\m{o,T;H)nc°([o,T}-y) + ||^r||//i(o,T;H r )nc ([o,T];Vr) < C2 u r )\\n ■ (2-15) 
Next, a comparison argument in (12.91) shows that also 

\\Aw\\ L 2 {0jT . H) < C 3 \\(u,u r )\\ H . (2.16) 
Now we invoke [1, Theorem 3.1, p. 1.79] with the specifications 

A = -A, g = yir, V = 2, r = 0, s = 3/2, 
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to conclude that 



o 



whence it follows that 



w(t)\\ 2 H3/2(a) dt < C 4 [ T (\\Aw(t)\\ 2 H + \\w r (t)\\ 2 Vr ) dt, (2.17) 



\\w\\i*(p,T;H*/*{n)) < C 5 \\{u,u r )\\n- (2.18) 
Hence, by the trace theorem, we have that 

\\d a w\\ L 2 {0>T . Hr) < Ce\\{u,u r )\\u, (2.19) 

so that, by comparison in the equation resulting from (12.101) . we obtain 

||A r w;r||L2 (E) < \\(u,u r )\\ H (2.20) 

and consequently 

\\wr\\L2(0,T;H2(T)) < Cg \\(u, U r )\\ H . (2.21) 

Now, owing to standard elliptic estimates, we infer 

IHU 2 (o,tw 2 (^)) - c 7\\(u,u r )\\ n . (2.22) 
This concludes the proof of the assertion. B 

Remark 3: It follows from (ii) in Theorem 2.2 that for zero initial data the solution 
operator (u, ur) (w, wr) is a continuous linear mapping from Ti into y . 

While it cannot be expected that the solution to the linear system (I2.9l) - (l2.1ip is bounded, 
we now establish a boundedness result for the solution to the nonlinear state system (11.21) - 
(11.41) that will be of key importance in the subsequent analysis. To this end, we need the 
following assumption: 

(A5) y e L°°{Q) , y 0r G L°°(r) , and it holds 

< essinf yo(x), ess sup yo(x) < 1, 

< essinf y 0r (x), esssup yo r (x) < 1. (2.23) 



Lemma 2.3 Suppose that the assumptions (A2)-(A5) are satisfied. Then there are 
constants < r* < r* < 1 , which only depend on Q , T , y , y 0r , f , g , and R , such 
that we have: whenever (y,yr) = S(u,ur) for some (u,ur) G U then it holds 

r* < y < r* a. e. in Q, r* < yr < r* a. e. in E. (2.24) 

Proof: Let (u, ur) G U be arbitrary and (y, y-p) = S(u, u^) . Then we have 

IM|l°°(q) + ||wr||L°°(E) < R- 
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By virtue of ( 12. 2 j) and (I2.23j) . there are constants < r* < r* < 1 such that 



r* < min < essinf yo(x) , essinf yo F (x) > , (2.25) 

r* > max < ess sup yo(x) , ess sup yo T {x) \ , (2.26) 

max{/'(r) + R, g'(r) + R} < VrG(0,r,), (2.27) 
min{/'(r) -R, g'(r) - R} > Vre(r*,l). (2.28) 

Now define w := (y — r*) + . Clearly, we have w G V and w\r G Vp. We put Wr '■= w\r 
and test (jl.2p by w . Thanks to (I2.26P , we readily see that 

= \\\w{T)\\ 2 H + £\\Vw(t)f H dt 

+ 1\\MT)\\h t + / HVr^r(t)|lH r dt + $, 
where, owing to (12T27D and (l2T28|) . 

* := / [ (f(y)-u)wdxdt + / [{g'{y r )-Ur)w r dTdt > 0. 
jo Jo ir 

In conclusion, w = (y — r*) + = 0, i.e., y < r* , almost everywhere in Q and on S. The 
remaining inequalities follow similarly by testing ( II. 2p with w := — (y — r*) _ . B 



Observe that in view of (A2) and of Lemma 2.3, we may (by possibly choosing a larger 
-frj* ) assume that also 

max {max {||/%)IU~ (Q) , ||<?%r)IU~( S )}} < ^ , (2.29) 
whenever (y,yr) = S(u,Ur) for some (u,ur) G U . 

Remark 4: Lemma 2.3 entails that the singular components in the state equations of 
the control problem (CP) are only active in a domain of arguments where they behave 
like standard bounded smooth nonlinearities. As a consequence, we could use classical 
differentiability results to see that both / and g generate three times continuously dif- 
ferentiable Nemytskii operators on suitable subsets of L°°(Q) and L°°(E), respectively. 
From this point, it would in principle be possible to derive the subsequent differentia- 
bility results for the control-to-state mapping by using the implicit function theorem. A 
corresponding approach was taken in Chapter 5 in [16] for the case of standard Neumann 
boundary conditions not involving dynamic terms or the surface Laplacian. Here, we 
prefer a direct approach which, while being slighthly longer and possibly less elegant than 
the use of the implicit function theorem, has the advantage of yielding the explicit form 
of the corresponding derivatives directly. 
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With Lemma 2.3 at hand, we are now able to improve the stability estimate ( 12. 8 j) from 
Theorem 2.1. 

Lemma 2.4 Suppose that the general assumptions (A2)-(A5) are satisfied. Then there 
is a constant K% > , which only depends on Q , T , f , g , and R , such that the 
following holds: whenever (u\, u± r ), (u±, U2 r ) E U are given and (t/i, 2/i r ), (2/2, J/2 r ) E y 
are the associated solutions to the state system f 1 1 . 2 j) — ( TTT4|) . then we have 

WiVuVlv) ~ (Z/2,3/2r)||;y < Kl ||(«i,Wi r ) - (« 2 ,«2r)ll«- ( 2 - 30 ) 

Proof: In the following, we denote by Cj , i E IN , positive constants that only depend 
on the quantities mentioned in the assertion. We subtract the state equations (jl.2p - (jl.4p 
corresponding to ((ttj, u ir ), (yi, y ir )) , i = 1,2, from each other and multiply the equation 
resulting from flT5]) by - y 2 ) • Putting u = u x -u 2 , u T = u lv - u 2r , y = ?/i - y 2 , 
and yr = 2/i r — V2 r , we have for all t E [0, T] the estimate 

/7 y t 2 dxd S + lf \Vy(t)\ 2 dx + ft |% r | 2 drds + ~ f |V r y r (t)| 2 dr 

[ [\f'(yi)-f(y2)\\yt\dxds+ [ [ \g'(y lT ) - g'(y 2T )\\d t y r \dTds 
Jo Jn Jo Jr 



< 



+ / /"|u| |j/t|drds + / / |u r | \d t yr\ dTds. (2.31) 
jo Jn Jo Jr 

Now observe that Lemma 2.3 (see also (12.291) ) and the mean value theorem imply that 
there is some constant C\ > such that 

l/'G/i) - /' (2/2)| < Ci|y| a.e. in Q, |p'(yi r ) - ^(|/2 r )| < Ci \yr\ a.e. in S. 
Hence, it follows from Young's inequality and ( 12 .8p that 

\\y\\m(o,T;H)c\c°{[o,T}-y) + \\yr\\m{o,T;H T )c\c°{[o,T}-y T ) < C2 ||(w,ur)||w ■ (2.32) 

From this point we may continue as in the proof of Theorem 2.2 after proving the estimate 
(I2.15P : indeed, by the arguments used there, we can repeat the estimates ( 12.161) to (I2.2ip 
with (w,wr) replaced by (y,yr) to come to the conclusion that (I2.30p is satisfied. This 
concludes the proof of the assertion. B 

3 The optimal control problem 

We now consider the optimal control problem (CP) under the general assumptions (Al)- 
(A4). 



3.1 Existence 

We have the following existence result. 
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Theorem 3.1 Suppose that the general assumptions (A1)-(A4) are fulfilled. Then the 
optimal control problem (CP) admits a solution. 

Proof: Let {(u n ,u nr )} C be a minimizing sequence for (CP), and let (y n ,y nr ) = 
S(u n ,u nr ) , n G IN. By virtue of the global estimates (12. 6 p and (I2.24p . we may assume 
(by possibly selecting a suitable subsequence again indexed by n ) that there are functions 
(u,u~r) £ X and (y,yr) £ y, such that 

u n — > u weakly-* in L°°(Q), 

u nv — > ur weakly-* in L°°(£), 

y n ^y weakly-* in H 1 ^, T; H) PI L°°(0, T; V) fl L 2 (0, T; H 2 (£l)) H L°°(Q), 

yn T ^yv weakly-* in if 1 (0, T; if r ) H L°°(0, T; Vr) H £ 2 (0, T; if 2 (r)) fl L°°(S) . 

In particular, we have 

d n y n -> d n y, A r y nr -> A r y r , both weakly in L 2 (S). 

Clearly, («, w-) G i/ a d • Moreover, we infer from standard compact embedding results 
(cf. [151 Sect. 8, Cor. 4] that, in particular, 

y n -> y strongly in C°([0,T]; H), (3.1) 
y nr -> y r strongly in C°([0, T]; # r ). (3.2) 

But then we can conclude from the Lipschitz continuity of f 2 and g' 2 (see (A2)) that also 

&(Vn) -> strongly in C°([0, T];H), 

fl^(j/n r ) fi^r) strongly in C°([0, T]; #-), 

while (12. 7p and (A2) allow us to deduce that 

f[(yn)^C weakly in L 2 (0,T;tf), 
9i(yn r ) -> £r weakly in L 2 (0, T; if r ). 

for some weak limits £ and Since /i and are convex (so that f[ and are 
increasing), the weak convergences above, along with (I3.1I) - (I3.2I) . imply that £ = f[{y) 
a.e. in Q and £r = fi'i(yr) a-e. on S, due to the maximal monotonicity of the operators 
induced by f[ on L 2 (Q) and g[ on L 2 (£) (see, e.g., [H, Prop. 2.5, p. 27]). At this 
point, we may pass to the limit as n — > oo in the state system (jl.2p ~ p.4p (written for 
(y n ,y nr )i K)V)i n e ^° conclude that (y,y~r) = «£(«,«-), that is, the pair ((w, Wr), 
(y,y-p)) is admissible for (CP). It then follows from the lower sequential semicontinuity 
of the cost functional J that (u, Mr) is in fact an optimal control for (CP). - 

3.2 Differentiability of the control-to-state operator 

Suppose now that (u, ur) £ is a local minimizer for (CP), and let (y, y r ) = S(u, Ur) 
be the associated state. We consider for fixed (h, h r ) e X the linearized system: 
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£t-A£ + f"(y)Z = h a.e. in Q, (3.3) 

<9 n £ + dtfr - A r £r + g"(Vr) £r = h r , £ r = £|r, a. e. on £, (3.4) 

f(-,0) = a.e. in fi, fr(-,0) = a.e. on T. (3.5) 



By Theorem 2.2 the system fl3.3l) - fl3.5l) admits for every (h, hr) G H (and thus, a fortiori, 
for every (h,hr) E X) a, unique solution (£,£r) G 3^, and the linear mapping (h,hr) i-> 
(£, £r) is continuous from "H into 3^ and thus also from X into 3^- 
We have the following differentiability result. 

Theorem 3.2 Suppose that the assumptions (A2)-(A5) are satisfied. Then we have the 
following results: 

(i) Let (u,ur) G U be arbitrary. Then the control-to-state mapping S, viewed as a 
mapping from X into y , is Frechet differentiable at (u, u~r) , and the Frechet derivative 
DS(u,ur) is given by DS(u,u r )(h, hr) = (£, £r)> where for any given (h,hr) G X the 
pair (£, £r) denotes the solution to the linearized system f)3.3l) - fl3.5l) . 

(ii) The mapping DS : U — > C(X,y), (u, Up) H- DS(u,ur) is Lipschitz continuous on 
U in the following sense: there is a constant K\ > such that for all (u,ur), (u,Tir) G U 
and all (h, hr) G X it holds 

\\{DS(u,u r ) - DS(u,u r )){h,h r )\\y < K\ \\ (u, u r ) - (u, u r ) \\ H II (h, h r ) \\ H . (3.6) 

Proof: We first show (i). To this end, let (u,ur) G U be arbitrarily chosen, and let 
(y,yr) = <S(u,ur) G 3^ denote the associated solution to the state system. Since U is an 
open subset of X , there is some A > such that for any (h, h r ) G X with || (h, hr) \\x < \ 
it holds (u,u~r) + {h, h r ) G U . In the following, we consider such variations (h, h r ) G X , 
and we denote by (y , yr) the solution to the state system fll.2p - fll.4l) associated with 
(u, ur) + (h,hr). Moreover, we denote by (^ ft ,^r) the unique solution to the linearized 
system fl3.3p - fl3.5l) associated with (h, hr) . We also denote by Ci , i G IN , positive 
constants that depend neither on the choice of (h, hr) G X with \\(h, hr)\\x < A nor on 
t G [0, T] . 
Now let 

v h :=y h -y-C h , t£ := j£ - - & 

Obviously, we have (v h ,Vr) G y. Since the linear mapping (h,hr) ((, h ,(,r) * s by 
Theorem 2.2 (ii) continuous from X into y , it obviously suffices to show that there is an 
increasing function G : [0, A] — > [0, +oo) which satisfies lim r \^ G(r)/r = and 



\\(v\v£)\\ y < G(||(Mr)||«) • (3-7) 
Apparently, v h is a solution to the initial-boundary value problem 

d n v h + d t v h T - A r v? + </'(^) " g'ivv) ~ g"(Vr) Cv = a. e. on E, (3.9) 
u h (-,0) = a.e. in ft ^(-,0) = a.e. on T. (3.10) 



12 



Optimal control of an Allen-Cahn equation 



Next, we observe that it follows from Taylor's theorem and from the global estimate (I2.24p 
(cf. also ( I2.29P ) that almost everywhere on Q we have 

f'(y h ) - f(V) - f"(y) e = f"(V) v h + $ fc , (s.n) 

with some function $ fe G L°°(Q) such that, almost everywhere in Q, 

\$ h \<l max \f^( r )\\y h -y\ 2 <^f\y h -y\\ (3.12) 
By the same token, there is a function $r £ L°°(S) such that, almost everywhere on S, 

g'ivr) - g'ivr) - g"(yr) & = g"(yr) wf + (3.13) 

where 

|*?| < ^bNyrl 2 . (3.14) 

Hence, with Cl := /"(£) 6 c 2 := g"(y r ) G £°°(£), u := G L 2 (Q), and 

u r := — $r G £ 2 (X) , we see that the system fl3.8IHI3.10l) satisfied by (v h ,v^) has exactly 
the structure of the system fl2.9l) - fl2.1ip . It therefore follows from (I2.12j) in Theorem 2.2 
that 

Mv h ,v$)\\y < C^K&^Wh. (3.15) 

Now observe that, owing to the embedding V C L 4 (Q) and to (I2.30P in Lemma 2.4, we 
have 

||*iia W) < C 2 f T f \y h -y\ 4 dxdt = C 2 [ T \\y h (t)-y(t)\\Un)dt 
Jo Jn Jo 

< C 2 T \\y h - y\\* c0mnv) < C 3 \\(h,h r )\\ 4 n . (3.16) 
Similar reasoning shows that also 

ll^ll!^) < C A \\{hM)t H - 

In conclusion, (13. 7p is satisfied for the function G{r) = Ci(y/Qi + \[Ca) r<1 , which con- 
cludes the proof of assertion (i). 

Next, we show the validity of (ii). To this end, let (u, Wp) G U be arbitrary, and let 
(k, kr) G X be such that (u + k,u r + k r ) G U . We denote (y k , y*) = S(u + k,u r + k r ) 
and (y,y~r) = S(u,ur), and we assume that any (h,hr) G X with \\(h, hr)\\x = 1 is 
given. It then suffices to show that there is some L > , independent of (h, h?) , (u, «r) 
and (k, kr) , such that 

ll(£*,#)-(£,fr)b < L\\(k,k r )\\ H , (3.17) 

where (£ fc , £*) = DS(u + k, u r + kr)(h, hr) and (£, £r) = DS(u, ur)(h, hr) ■ For this 
purpose, we denote in the following by , z G IN, positive constants that neither depend 
on (u, u r ) , (k, kr) , nor on the special choice of (h, h r ) G X with \\(h, h r )\\x = 1 . 
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To begin with, observe that from part (i) it easily follows that (w,Wr) '■= (£ fc ,£r) — 
(£, £r) G y solves the initial-boundary value problem 

w t -Aw + f"(y) w = - e (f"(y k ) - f\y)) a. e. in Q, (3.18) 
d n w + d t w r - A r w r + g"(y r ) w T = - (g"(y k ) - g"(y r )) a. e. on £, (3.19) 
io(-,0) = a.e. in SI, i« r (-,0) = a.e. on T. (3.20) 

Hence, it follows from Theorem 2.2 that 

\\(w,w r )\\ y < k\ (\\e(f"(y k )-f"m\\L>( Q ) + UHg"(yr)-g"(m\L^)) • (3.21) 

Now, by Taylor's theorem and (12.291) . we have almost everywhere in Q (on E, respec- 
tively) 

l/V) - f"(y)\ < K{ \y k - y\ and \g"{y k ) - g"{y T )\ < K{ \y k T - y r \ . (3.22) 

At this point, we recall that IA is obviously a bounded subset of "H . Since (u + k, u + k-p) G 
U and || (h, hr) \\x = 1 , we thus can infer from (12. 29ft and from the estimate (12. 12ft in 
Theorem 2.2 that (£ fc ,£p) is bounded in y independently of (k, k?) , (u,ur), and the 
choice of (h,hr) G X with \\(h, hr)\\x — 1- Using the embedding V C L 4 (f2) and 
Lemma 2.4, we therefore have 

U k (f"(y k )-f"(y))\\l HQ) < K 2 f f \e\ 2 \y k -y\ 2 dxdt 

Jo Jn 

< K 2 [ T Wemhw \\y k (t)-y(t)\\% in) dt 

Jo 

< K 3 \\(y k -y,y k -y T )\\ 2 y < K A \\ (k, k T ) \\ 2 n . (3.23) 

Since an analogous estimate holds for the second summand in the bracket on the right- 
hand side of (I3.2ip . the assertion follows. a 

Remark 5: Notice that we could not establish Frechet differentiability of S on T-L; 
in fact, we only were able to show differentiability on the bounded subset U of X . 
In particular, the boundedness of (u, ur) in X was an indispensable prerequisite for 
proving Lemma 2.3 and the global estimate (1 2 . 2 9 j) . which in turn was fundamental for 
the derivation of the differentiability result. This will below lead to a so-called two-norm 
discrepancy in the derivation of second-order sufficient optimality conditions, i.e., we will 
have to work with two different norms. 

With Theorem 3.2 at hand, it is now straightforward to derive the standard variational 
inequality that optimal controls must satisfy: indeed, it follows from the quadratic form 
of J and from the chain rule that the reduced cost functional 

J{u, u r ) := J(S(u, u T ), (u, u T )) (3.24) 

is Frechet different iable at every (u, Ur) G U with the Frechet derivative 

DJ(u,u r ) = D {ytyr) J(S(u,u T ), (u,u r ))oDS(u,u r ) + D {UjUr) J(S(u,u r ), (u,u r )) , (3.25) 
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and, owing to the convexity of W a d , we have for every minimizer (u, Mr) G U&d of J in 
W a d that 

DJ{u,u T ){v-u,v T -u T ) > V(t),!j r )6^. (3.26) 

Identification of the expressions in (I3.25P from (11. ip and Theorem 3.2 yields the following 
result: 

Corollary 3.3 Let the assumptions (A1)-(A5) be satisfied, and let (u, Ur) G W a d be an 
optimal control for the control problem (CP) with associated state (y, yr) = S(u, ur) G y . 
Then we have for every (v, Vr) G £4d 

hi l{y-z Q )£dxdt + fol /(yr-^Krdrdt 

Jo Jn Jo Jr 



+ /3 3 / (y(-,T)-z T )a-,T)dx + p 4 / (y r (-,T)-zr, T )^(-,T)dT 
Jn Jr 

+05 I fu{v-u)dxdt + p 6 I fu T {v v -u T )dYdt > 0, (3.27) 
Jo Jn Jo Jr 

where (£,£r) G y is the unique solution to the linearized system f l3.3l) - f)3.5p associated 

with (h, hr) — (v — u, vr — Mr) ■ 

3.3 First-order necessary optimality conditions 

We are now in the position to derive the first-order necessary optimality conditions for 
the control problem (CP). 

Theorem 3.4 Let the assumptions (A1)-(A5) be satisfied, and let (w, «r) G be an 
optimal control for the control problem (CP) with associated state (y,yr) = S(u,u-p) G y . 
Then the adjoint state system 

- Pt -Ap + f"{y) p = Pi(y- zq) a. e. in Q, (3.28) 
d n P ~ d t pr ~ A r pr + g"{yr) Pr = Pi (Vr - Zn) a. e. on £, (3.29) 
P( • , T ) = Pz(V{ ■ , T) - z T ) a. e. in ft, 

pr{ ■ , T) = /3 4 (yr{ ■ , T) - z r>T ) a. e. on T, (3.30) 
has a unique solution (p,pr) G y , and for every (v,Vr) G W a d we have 

(p + P 5 u)(v -u)dxdt + / (p r + p 6 u r )(v r -u r )dTdt > 0. (3.31) 




o Jn Jo Jr 



Proof: First observe that the system (13. 280 - 03.300 is a linear backward-in-time parabolic 
initial-boundary value problem, which after the time transformation t i— > T — t takes the 
form of the system fl2.9p - fl2.lip provided we put 

d{x,t) : = f"(y(x,T-t)), c 2 {x,t) := g"(y r (x,T - t)), 

u(x,t) := fix (y - z Q )(x,T- t), u T (x,t) := P 2 (Vr - z^)(x,T - t), 

w (x) := p 3 (y(x,T) - z r (x)), w 0r (x) = p 4 (y r (x,T) - z v>T {x)) . 
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Obviously, ci G L°°{Q), c 2 G L°°(E) , u G £ 2 (Q), «r G £ 2 (£), u> G i? 1 ^) , and 
u>o r ^ • Hence, by virtue of Theorem 2.2, the transformed system admits a unique 

solution (w,Wr) G 3^, so that (p,pr)(^,^) : = (w,Wr)(x, T — t) is the unique solution to 
the adjoint system, and (p,pr) G y . 

At this point, we may perform the standard calculation, using repeated integration by 
parts and the systems fl3.3p -f l3"3j) and fl3.28p - fl3.30p . which shows that 

Pi[ f{y-z Q )£dxdt + fof [{yr-z^frdTdt 
Jo Jn Jo Jr 

+ (3 3 [ (y{-,T)-z T )C{-,T)dx + (3 4 [ {y r ( ■ , T) - z r , T ) M • , T) dT 
Jo, Jr 

phdxdt+ [ [p r h r dTdt, (3.32) 




'0 Jn Jo Jr 

so that f EPTj) follows from fl3T27j) . 



Remark 6: (i) It follows from the above considerations that the Frechet derivative 
DJ{u,ut) G C(X,y) can be identified with the pair (p + f3 5 u,pr + P&uy) in the sense 
that, with the standard inner product (• , • )% in the Hilbert space "H, we have 

DJ(u,u r )(h,h r ) = ((p + f3 5 u,p r + f3 6 u r ),(h,h r )) H W(h,h r )eX. (3.33) 

(ii) If /3 5 > and (3 e > then it follows from standard arguments that the condition 
(13. 3 1 p can be given a pointwise interpretation in the following sense: we have almost 
everywhere in Q (on £ , respectively) that 

{u 2 (x,t) if u 2 (x,t) < -f3^p(x,t) 

-f3^ l p{x,t) if ui(x,t) < -I3^ 1 p(x,t) < u 2 (x,t) , 
Ui(x,t) if ui{x,t) > -P^p{x,t) 

(u 2r (x,t) if u 2r (x,t) < -(3Q l p r (x,t) 

-(3Q l p T (x,t) if u lr {x,t) < -Pq 1 p r {x,t) <u 2r {x,t) , (3.34) 
u lr (x,t) if u lr {x,t) > -Pq 1 p r (x,t) 

where Ui, u 2 , U\ r , u 2r represent the control constraints defined in (Al). 



3.4 The second-order Frechet derivative of the control-to-state 
operator 

Since the control problem (CP) is nonconvex, the first-order necessary optimality con- 
ditions established in the previous section are not sufficient. However, it is of utmost 
importance, e.g., for the numerical solution of (CP), to derive sufficient optimality con- 
ditions. For this purpose, it is necessary to show that the control-to-state mapping is 
twice continuously Frechet differentiable. We have the following result. 



16 



Optimal control of an Allen-Cahn equation 



Theorem 3.5 Assume that in addition to the general conditions (A1)-(A5) we have: 

(A6) /,<?eC 4 (0,l). 

Then we have the following results: 

(i) The control-to- state operator S is at any (u, ur) E U twice Frechet differentiable, 
and the second Frechet derivative D 2 S(u,ur) E £(X,£(X,y)) is defined as follows: if 
(h, hr), (k, kr) E X are arbitrary, then D 2 S(u,ur)[(h, hr) , (k,kr)] —: (t),t)t) E y is the 
unique solution to the initial-boundary value problem 

Vt-A V + f"(y) V = -f® (y) ^ a . e in Q, (3.35) 

d n r] + d t r]r - A r r/r + g"(yr) Vr = -0 (3) (Vr) <Pr a. e. on E, (3.36) 

r/ ( • , 0) = a. e. in Q, r/r( ■ , 0) = a. e. on T, (3.37) 

where we have put 

(y,y r ) = S(u,u r ), (0,0 r ) = DS(u,u T )(h,h T ), (ifi,ipr) = DS(u,u r ){k,k r ) . (3.38) 

(ii) The mapping D 2 S : U —> £(X, £{X, y)) , (u,ur) D 2 S(u,ur) , is Lipschitz con- 
tinuous on U in the following sense: there exists a constant > such that for every 
(u, u r ), (u, Ur) E U and all (h, hr), (k, k r ) E X it holds 

\\(D 2 S(u, u r ) - D 2 S(u, u T ))[(h, h r ), (k, k r )}\\y 

< K* 5 II (u, u r ) - («, u T ) \\ H II (h, h T ) \\u II (k, k r ) \\ n ■ (3.39) 

Proof: We first prove part (i) of the assertion. To this end, let (u, u r ) E U be fixed, 
(h, hr), (k, kr) E X be arbitrary, and (y, y~r), (<p, 4>r), ("0> V'r) E y be defined as in (13.381) . 
Then, with 

Cl :=f"(y)EL°°(Q), c 2 := g"(y~r) E L°°(E), u := -/W(y) 0^ E L 2 (Q), 
«r:=-S (3) (yr)0r^reL 2 (E), 

the system (I3.35j) - (j3.37j) takes the form (j2.9B - fl2.llj) and thus enjoys a unique solution 
(r),r]r) E y. Moreover, by fl2.12j) we have the estimate 

\\{v,Vr)\\y < C(\\f®(y)M\\v {Q) + lk (3) (yr)0r^r||L 2( E)) • (3.40) 

In the remainder of the proof of part (i), we denote by Cj , i E IN, positive constants that 
do not depend on the quantities (h,hr), (k,kr), and (u, Ur) ■ Using (I2.29|) and (12.301) . 
and invoking the embedding V C L 4 (f2) , we find that 

\\f {3) (y)<i>Hl HQ) < Kfff |0| 2 H 2 d*dt < d f \\m\\i^)\\m\\i^t 

Jo Jn Jo 

< C2 U\\ 2 c 0{[ o,t];V) UWl^nv) < C 3 \\(h,hr)\\ 2 n \\(k,kr)f n , ( 3 - 41 ) 

where the validity of the last inequality can be seen as follows: by definition (recall (I3.38P ) 
((f), (fir) is the unique solution to the linear problem (I3.3j) - (j3.5p with zero initial conditions, 
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and thus we can infer from Theorem 2.2 (see, in particular, (I2.12p ) that ||(0, 0r)||y < 
C \\(h,hr)\\u- By the same token, we conclude that || (-0, Vr) ||y < C \\ (k, kr) \\n • The 
asserted inequality therefore follows from the definition of the norm of the space y , and 
we obtain from similar reasoning that also 

lb (3) (yr)0r^r|U 2(E) < C 4 \\{h, h r )\\ H \\(k,k r )\\ H - 

In particular, it follows that the bilinear mapping X x X y , [(k, kr), (h, hr)] !->■ (rj, rjr) , 
is continuous. 

Now we prove the assertions concerning existence and form of the second Frechet deriva- 
tive. Since IA is open, there is some A > such that (u + k,Ur + kr) G U whenever 
\\{k, kr)\\x < A. In the following, we only consider such perturbations (k,kr) EX. Then 
for (y k ,y£) = S(u + k,u r + k T ) the global estimates flZD, fl2T29|) . flODD are satisfied. 
Without loss of generality, we may also assume that 

max {||/ (4 V)|U~(0), \\9 {A) (y k r)\\L^)} < K{ whenever ||(Mr)|U < A. (3.42) 

After these preparations, we observe that it suffices to show that 

LD«S(m + k, ur + kr) — DS(u, ur) — D 2 S(u, Ur)(k, kr 



c{x,y) 

\(h h„\\ (h h~\\\ 

< G(\\{k,kr)\\u) , (3-43) 



sup || (DS(u + k, ur + kr) — DS(u, ur) — D 2 S(u, ur)(k, k r )) (h, hr) 
\\(h,h r )\\ x =i 



with an increasing function G : (0, A] — > (0, +00) that satisfies lim r ^o G(r)/r = 0. 
To this end, let (h,hr) G X be arbitrary with ||/i||i°°(Q) + ||/ir|U°°(E) = 1- We put 
(p, pr) = DS{u + k,yr + kr)(h, hr) , define the pairs (0, 4>r), ^r) G y as in f)3.38p . and 
put 

(w, w r ) := (p - - 77, p r - 0r - Vr)- 
Then, according to f)3.43p . we need to show that 

\\(w,wr)\\y < G(\\(k,kr)\\ n ) . (3.44) 

Now, invoking the explicit expressions for the quantities defined above, it is easily seen 
that (w, Wr) is a solution to the linear initial-boundary value problem 

w t — Aw + f"(y) w = a a. e. in Q, (3.45) 
d n w + d t wr — ArWr + g" (yr) wr = &r a. e. in E, (3.46) 
w{-,0) = a.e. in Q, w r (-,0)=0 a. e. on T, (3.47) 

where we have put 

o- := -p (/" (y k ) - f" (y)) + / (3) 0/)<^, 

*r := -Pr (sfW - 9"(Vr)) + <? (3 %r) </>r>r • (3.48) 
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In view of (I2.29p . and since it is easily checked that (a, op) G "H, we may again invoke 
the estimate (I2.12p in Theorem 2.2 to conclude that (I3.44p is satisfied if only 

\\(cr,a r )\\ n < G (\\(k,k r )\\ H ) . (3.49) 

Applying Taylor's theorem to and recalling (I3.38p . we readily see that there is a 
function Uf G L°°(Q) such that 

f"(y k ) - f"(y) = f m (y) (y k -y-^) + f (3) (y) 1> + w/ (y k - y? a. e. in Q . (3.50) 

Hence, we have that 

o =- P f {3) (y)(y k -y-^) - ^f (3) (y) (p - <f>) - P oj f (y k -y) 2 . (3.51) 

Now recall that from the proofs of Frechet differentiablity and of the Lipschitz continuity 
of the Frechet derivative (see, in particular, the estimates fl3.6l) - f)3.16l) and fl3.17l) - fl3.23l) . 
respectively) it follows that 

\\{y k -y-^^-y v -^)\\ y < & ||(Mr)||L 

||(p-0,pr-0r)b < C 2 \\(k,k r )\\ H . (3.52) 

Moreover, we can infer from Lemma 2.4 that 

\\(y k -y,yr-yr)\\y < K*\\(k,k r )\\ n , (3.53) 

and it follows from Theorem 2.2 that p is bounded in y by a positive constant that is 
independent of (k, k r ), (h, h r ) G X with ||(fc, fcr)||;t < A and /ir)IU = 1 • 
Finally, we conclude from Remark 3 that with a suitable constant C3 > it holds 

||(V^r)b < C 3 || (Mr) Ik- (3.54) 

After these preparations, and invoking Holder's inequality and the embeddings V C L 4 (Q) 
and V C L 6 (f2) , we can estimate as follows: 

\W\\h iQ) < cJ T [ (|p|V-y-<A| 2 + H 2 Ip-0I 2 + \p\ 2 \y k -y\ 4 ) dxdt 

Jo Jn 

< ^ T (llK*)lli 4( n)ll(/-^-V')(t)ll!^) + IW)lli* ( n)llp(*)-^)lli*(n))d* 

T 



+ C4 / \\ P (t)\\i 6in) mt) -y(t)\\h(n } dt 
Jo 

< c 5 max (\\p(t)\\v\\(y k -y-m)\\v + \mt)\\ 2 v\\p(t)-m\ 2 v 

+ \\p(t)\\ 2 v\\y k (t)-mrv) 

< C 6 \\(k,k r )\\ 4 n . (3.55) 

By the same reasoning, a similar estimate can be derived for ||cx r || L 2( S ) , which concludes 
the proof of the assertion (i). 
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Next, we prove the assertion (ii). To this end, suppose that (u, Ur) E U and that (h, hr) 
and (k, kr) are arbitrarily chosen in X , and let (5, 5r) E X be arbitrary with (u + 5, ur + 
$r) E X . In the following, we will denote by Cj, j 6 IN, positive constants that do not 
depend on any of these quantities. We put 

(y 5 ,Vr) = S(u + 5,u r + 5r), (y,yr) = S(u,u r ), ((j>, 0r) = DS(u,u T )(h, h T ), 
(V, Vr) = DS(u, u v )(k, k r ), (4, 4r) = DS(u + 5,u r + 5 r )(h, h r ), 
(^, 4) = DS(u + 5,u r + 5 T )(k, kr), (v, Vr) = D 2 S(u, u T ) [(h, h v ), (k, kr)], 
(4, 4) = D 2 S(u + 5,u r + 5 r ) [(h, h v ), (k, k r )\ . 

From the previous results, in particular, (12.301) and (13.61) . we can infer that there is a 
constant C\ > such that 

\\(<f>,<f>r)\\y+\\(4,4r)\\y < Ci||(Mr)||«, 
ll(^r)b+||(^#)b < C x \\ (k, k r ) \\ H , 
\\(v,Vr)\\y + \\(v S ,Vr)\\y < ||(Mr)IMIto 
\\(y 5 ,y 5 r)-(y,yr)\\y < CilKWUw, 
\\(4,4r)-(<PAr)\\y < 0,11(5, 5r)\\n\\(h,hr)\\n, 
\m 5 ,4)-(^^r)\\y < 0,11(5, 5r)\\n\\(k,kr)\\n. (3.56) 

Now observe that (w,wr) = (4,4) ~~ {VjVv) satisfies the linear initial-boundary value 
problem of the type fl231)- fl2TTB 

Wt — Aw + f"(y) w = a a. e. in Q, (3.57) 
d n w + d t wr — A r wr + g"(yr) wr = &r a - e. on S, (3.58) 
w( ■ , 0) = a. e. in Q, wr(-,0) = a. e. on T, (3.59) 

where we have put 

° = -v 5 (f"(y s ) - f'm - (/ (3 V) 4^ - f {3) (y)<Pi>) , 

^ = -4(g"(y 5 r) - g"(yr)) - (g {3) (y S r ) 4r4 ~ 9 {3) (yr) 0r>r) • (3.60) 
From Theorem 2.2 it follows that 

\\(w,w r )\\y < d\\(a,a r )\\ n , (3.61) 
so that it remains to show an estimate of the form 

\\(a,a r )\\n < C 2 \\(5, 5 r )\\ n \\(h, h r )\\ H \\(k, k r )\\ n ■ (3.62) 
Moreover, we can infer from (I2.30p and (13.61) that, almost everywhere in Q, 

W\ < K{ (\4\ \y s -y\ + \4\ U' 5 \ \y 5 ~y\ + \4\ \i> 5 - 1>\ + M \4 ~ <P\) ■ (3-63) 
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Hence, by (I3.56p . and using Holder's inequality and the embedding V C L 4 (f2) , 

[ T [ \v S \ 2 \y S -y\ 2 dxdt < [ T \\v S (t)\\l Hn) \\(y s -y)(t)\\LHn)dt 
Jo Jn Jo 

< C3||^ir c o ([ o !T];y) ||2/ 5 -y|r c o ([ o !T];y) < C 4 || (5,5 r )\\ 2 H \\(h, h r )\\ 2 H \\ (k, k r )\\ 2 H . 

(3.64) 

Similar reasoning yields 

U\^-n\h [Q) + \m<P 5 - <P)\\l HQ) < a||(5,5r)ll«ll(^^r)ll«ll(^MII«- (3-65) 

Moreover, we invoke f!3.56p . Holder's inequality, and the embeddings V C L 4 (f2) and 
H 2 (tt) C L°°(tt) , to conclude that 



T 




10? |VT 1/ - V\ 2 dxdt < f T \\(y s - y)(t)\\Un) 11^(011^(0 W(t)\\l^ d* 
i Jo 



'o Jn Jo 

— ^6 H^llcOaO.T^V) l|V'' 5 |lc' ([0,T];V r ) 1 1 2/ ^ — y|li 2 (0,T;H 2 (O)) 

< C7||(5,MII«II(^^)II«II(^MIIh- (3-66) 

Finally, we can estimate ||or||L2( S ) deriving estimates similar to fl3.63p - fl3.66p . which 
proves the validity of the required estimate f)3.62p . With this, the assertion is completely 
proved. a 

3.5 Second-order sufficient optimality conditions 

With Theorem 3.5 at hand, the road is paved to derive sufficient conditions for optimality. 
But, because the control-to-state operator S is not Frechet differentiable on ~H, we are 
faced with the two-norm discrepancy, which makes it impossible to establish second-order 
sufficient optimality conditions by means of the same simple arguments as in the finite- 
dimensional case or, e. g., in the proof of Theorem 4.23 on page 231 in [16]. It will thus be 
necessary to tailor the conditions in such a way as to overcome the two-norm discrepancy. 
At the same time, for practical purposes the conditions should not be overly restrictive. 
For such an approach, we follow the lines of Chapter 5 in [16], here. Since many of the 
arguments developed here are rather similar to those employed in [12], we can afford to 
be sketchy and refer the reader to [TB] for full details. 

To begin with, the quadratic cost functional J, viewed as a mapping on y xU , is obviously 
twice continuously Frechet differentiable on y x U , and for any ((y, yr), (u, u-p)) G y x U 
and any ((v, v r ), (h, h r )), ((w, wr), (k, k r )) e y x X it holds 

D 2 J((y,yr), (u,u T ))[((v,v r ), (h, hr)), ((w,w r ), (Mr))] 

= Pif [ vwdxdt + (3 2 [ [v r w r dTdt + (3 3 [ v(T)w(T)dx 
Jo Jn Jo Jr Jn 

+ p 4 [ v r (T) w r (T) dT + fi s [ [ hkdxdt + (3 e [ [h r k r dTdt. (3.67) 
Jr Jo Jn Jo Jr 
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Hence, it follows from Theorem 3.5 and from the chain rule that the reduced cost func- 
tional J is twice continuously Frechet different iable on U . Now let (m, Mr) 6 W be fixed 
and (h, hr), (k, kr) G X be arbitrary. In accordance with our previous notation, we put 

(y,y r ) = S(u,u r ), (0,0 r ) = DS(u,u T ){h,h T ), (i^,ipr) = DS(u,u T ){k, k T ), 
(r],T]r) = D 2 S(u,u r )[(h, hr), (k, k r )} ■ 



Then a straightforward calculation resembling that carried out on page 241 in [16], using 
the chain rule as main tool, yields the equality 

D 2 J(u,u r )[(h, h r ), (k, k r )] = D( y>yr) J((y,y T ), (u,u r ))(rj,rj T ) 

+ D 2 J((y, y r ), (u, «r)) [((0, r ), (h, h r )) , ((^, #), (k, k r ))\ . (3.68) 

Now observe that the first summand of the right-hand side of (13.681) is equal to the 
expression 

0i / (V ~ z Q )r]dxdt + fa / (y T - zz) r]r dT dt 
Jo Jn Jo Jr 

+ fo[ (V(T) - zr) V (T) dx + fa f (y r (T) - z v>T ) r, r (T) dT (3.69) 
Jn Jr 

and that (77, r/r) solves a system of the form (I3.3p - (l3.5p . with h replaced by —f^(y)4>ip G 
L 2 (Q) and hr replaced by — g^\y ^^r^r G L 2 (Y>). Since the calculation leading to the 
identity f!3.32j) also works for right-hand sides in L 2 (Q) x L 2 (S) , we can infer that 

D (y,yr) J ((yiyr), (u,u r ))(v,Vr) 

= ~ f [pf {3 \y)4>ijdxdt - [ T [ P rg {3 \yr)4>r^rdTdt, (3.70) 
Jo Jn Jo Jr 

where (p,pr) G y is the adjoint state associated with ((y, yr), (u, ur)) . Summarizing, we 
have thus shown that it holds the representation formula 

D 2 J(u,u r )[{h,hr),{h,h r )} = f ! (fa ~pf {3) (y)) \<p\ 2 dxdt 

Jo Jn 

+ f [(02-prg {3) (yr))\M 2 dTdt + fa [ |0(T)| 2 dx + fa [ |0r(T)| 2 dr 
Jo Jr Jn Jr 

+ A>N£ a(Q ) + 0e \\hr\\l H x) ■ (3.71) 

Equality (13.711) gives rise to hope that, under appropriate conditions, _D 2 J7(m, Mr) might 
be a positive definite operator on a suitable subset of the space H . To formulate such 
a condition, assume that (u, u r ) G is a given control with associated state (y, yr) = 
S(u,ur) G y and adjoint state p G y satisfying (I3.28I) - (I3.30I) . We then introduce for 
fixed r > the set of strongly active constraints for (u, Mr) by 

A T (u,u r ) := {(x,t) E Q : \p(x,t) + fau(x,t)\ > r} 

U{(x,t)GS: \p T (x,t) + 6 u T (x,t)\ > r} . (3.72) 
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Apparently it follows from (I3.3ip that, depending on the signs of p(x,t) + (35u(x,t) and 
of pr(x,t) + (3QUr(x,t), the control values u(x,t) and ur(x,t) , respectively, attain one 
of the constraint values. We now define the r— critical cone C T (u, ur) to be the set of all 
(h, hr) E X such that 

= if (x,t) E A T (u,u r ) 
h{x,t) { >0 if u(x, t) = U\(x, t) and (x, t) E~ A T (u, Ur) , 
<0 if u(x,t) = u 2 (x,t) and (x, t) E" A T (u, u r ) 

= if (x,t) E A T (u,u r ) 
hr{x,t) ^ >0 if ur(x, t) = u\ T {x, t) and (x, t) E" A T {u, w r ) . (3.73) 
<0 if Mr(x, t) = U2 r (x, t) and (x, t) E" A T (u, u r ) 

After these preparations, we can formulate the second-order sufficient optimality condition 
as follows: 



There exist constants 5 > and r > such that 

D 2 J[(u,u r )[{h,h r ),{h,h r )) > 5\\(h,h r )\\ 2 H V(h,h r )EC T (u,u r ), 

where D 2 J[{u ) ur) [(h, h r ), (h, h r )] is given by f)3.7ip with (y, y r ) = S(u, ur), 

((f), 4>r) = DS{u, ur)(h, h r ) and the associated adjoint state (p,pr) ■ (3.74) 

The following result resembles Theorem 5.17 in [T6] . 

Theorem 3.6 Suppose that the assumptions (A1)-(A6) are satisfied, and assume that 
the triple (u,ur) E (y,yr) = <S(u,ur) E y and (p,pr) E y fulfills the first- order 

necessary optimality conditions fl3.28l) - f)3.3ip . Moreover, assume that the condition f)3.74p 
is fulfilled. Then there are constants e > and o > such that 

J(u,u T ) > J(u,u r ) + a \\(u - u,u r - «r)||^ 

whenever (u,u r ) EU^ and \\(u,u r ) — (u,u r )\\ x < e. (3.75) 
In particular, (u, ■%) is locally optimal in the sense of X . 

Proof: The proof closely resembles that of Theorem 5.17 in [16], and therefore we can 
refer to [16]. We only indicate one argument that needs a bit more explanation. To this 
end, let (u, u r ) E be arbitrary. Since J is twice continuously Frechet differentiable 
in U , it follows from Taylor's theorem with integral remainder (see, e. g., Theorem 8.14.3 
on page 186 in |4]) that 

J(u,u r ) - J(u,u T ) = DJ(u,u r )(v,v r ) + - D 2 J(u, u r )[(v, v r ), (v, v T )} 

+ R J ((u,u r ),(u,u r )), (3.76) 
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with the remainder 

R J ((u,u T ), (u, u r )) 

= [ (l-s)(D 2 J(u + sv,u T + sv T )-D 2 J(u,u r ))[(v,v r ),(v,v r )}ds. 
Jo 

(3.77) 

A lengthy but straightforward calculation, based on the representation formulas (I3.67p - 
(I3.69P as well as on the Lipschitz estimates ( I2.30p . ( I3.6p . and ( I3.39p . reveals that 

\R J ((u,u T ),(u,u T ))\ < d / (l-s)s\\(v,v r )\\ 3 H ds 

Jo 

< C 2 \\(v,vr)\\ x \\(v,v r )\\ 2 H , (3.78) 

with global constants C\ > and C*2 > that do not depend on the choice of (u, u^) e 
W a( j. From this point, we can argue along exactly the same lines as on pages 292-294 in 
the proof of Theorem 5.17 in [16] to conclude the validity of the assertion. B 
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